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Abstract 

m ; 

The mass of the newly discovered pentaquark is calculated within the framework 
of the SU(3) Skyrme Model. Various estimates based on the model independent 
7-H ■ approach are compared with the model results and with the Chiral Quark Model. 

Our discussion shows that + is light with the mass of the order 1.5 GeV. 

m ■ 

O ' 1- Introduction 

^ : 

Ph! It seems that the new exotic spin 1/2 baryon of strangeness +1, + , is now well 

Q_i! established experimentally pQ , [2j , [2j , [I] . After the first report by Nakano, three other 

experiments have already confirmed the existence of + with mass Mq+ ~ 1540 MeV and 
very small width r e + < 25 MeV. It is interesting to note that such a state was predicted 
in the late 80-ies within the framework of the Skyrme Model. Indeed, Manohar ;5j and 
Chemtob in Ref.® mentioned that the SU(3) Skyrme Model 0,13,0 (SM) not only 
reproduced the spectrum of the lowest baryon multiplets 8 and 10, but also predicted the 
new, exotic states belonging to the higher SU(3) representations like 10 or 27. Such states 
were also seen in the KN scattering phase shifts calculated within SM in Refs. JT0] , \H\ ■ In 
1987 in Ref.[12j we presented the first estimate of the mass of the lightest member of 10 - 
today's + - with the result M®+ ~ 1530 MeV; in a striking agreement with the present 
experimental observations. Since the details of this estimate were never published, we 
think that it is worthwhile to recall shortly how this result was actually obtained. 

The real boost in this field was due to Ref . [THj published in 1997, where not only 
the mass of G + but also its width was calculated within the framework of the Chiral 
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Quark Model (xQM) [Tlj. In fact both SM and xQM have similar group theoretical 
structure and the main problem in both models consisted in estimating the average mass 
of the exotic antidecuplet, Mjq. The splittings within 10 are to a good approximation 
in the first order in the symmetry breaker, m s , identical as in the usual decuplet, i.e. 
proportional to the hypercharge Y, with proportionality coefficient being of the order of 
150 MeV. The problem of fixing Mj^ was solved in Ref.^3! by the assumption that the 
nucleon resonance N*(1710) belongs to 10. An immediate consequence of this assignment 
was that M e + ~ M^* — 150 = 1560 MeV. A more refined analysis of the antidecuplet 
splittings pushed the G + mass further down to the value of 1530 MeV [T3] . 

In our old estimate of the G + mass [T2] we used the second order perturbation theory 
in the strange quark mass m s . This allowed us to estimate the strange moment of inertia 
of the rotating soliton which is responsible for the 10 — 8 mass splitting. Therefore our 
analysis did not rely on any physical assumption assigning some existing nucleon resonance 
to 10. The price, however, was that the mass of + depended rather strongly on the value 
of m s and/or pion nucleon Y, n N term. 

2. Baryons in the chiral models 

Chiral models like xQM or SM are closely related to QCD [14] . One can formally 
integrate out gluons from the QCD Lagrangian, then the resulting nonlocal quark theory 
would necessarily respect chiral symmetry. One can approximate this complicated theory 
by a simple, chirally symmetric, quartic quark interactions. A model with such properties 
was formulated by Nambu and Jona-Lasinio however in a different context. NJL 
model exhibits an important phenomenon: spontaneous chiral symmetry breaking. As a 
result massless quark- ant iquark bound states emerge - pions, kaons and rj ((p a in short). 
Formally, one can now reexpress the interaction part of the NJL Lagrangian introducing 
auxiliary, composite meson fields. This is the Chiral Quark Model used in Ref.fTB*] and 
derived from QCD within the instanon model of the QCD vacuum [To] . 

Notice the absence of the explicit meson kinetic energy and meson self-interaction 
terms in xQM: they arise only from the quark loops when one integrates out the quark 
fields. The resulting Lagrangian can be organized in terms of a number of derivatives of 
the meson fields (gradient expansion) ^7j. A simple Ansatz for such a Lagrangian density 
was proposed by Skyrme 40 years ago JT5] : 

C= ^Tr(d^d^U) + ^Tr ([dJJ ,d v U rff) + r W z- (1) 

Here U = exp (i2tp a \ a /F n ), F^ = 186 MeV, e denotes the Skyrme parameter of the order 
4 — 5 and T W z stands for the Wess-Zumiono term [T5].|2*0]. 

Clearly, both xQM and SM are devised to describe meson physics at low energies. 
However, as realized by Skyrme [IB] and later by Witten j2Q] and collaborators pT] . there 
exists a solitonic solution - a nontrivial classical configuration of the meson fields, taken 
in the form of the hedgehog Ansatz: 

in-rP(r) n 

(2) 
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- which can be interpreted as a baryon. Further quantization of the rotational zero modes 
(both in the configuration space and in the flavor space), i.e. rotations of the meson fields 
U = AUqA^ parameterized by a time dependent SU(3) matrix A(t), provides quantum 
numbers corresponding to the different baryonic states [I|,[B],[H]- in the chiral symmetry 
limit the effective baryonic Hamiltonian takes the following form: 

#o = M cl + ±-S{S + 1) + i (c 2 (K) - S(S + 1) - ^pj . (3) 

Here S denotes baryon spin, C2(7Z) the Casimir operator for the SU(3) representation 
K = (p,q): 

C 2 {K) = l -{p 2 + q 2 +pq + 3(p + q)) . (4) 

The classical soliton mass M c \ [P] , and two moments of inertia [P] are functionals of 
the solitonic solution (J2J) and can be numerically calculated within the model considered. 
However, as already observed in Ref.[2I] in the SU(2) version of the Skyrme Model and 
further confirmed in the SU(3) case [22], the classical soliton mass is by far too large 
for the realistic set of parameters [22! . Therefore we adopt here a "model- independent" 
approach of Ref . j24j where M c \ and 2^2 are treated as free parameters. This approach may 
be justified by the large N c argument. Constants M c \, ii >2 ~ N c , hence the classical mass is 
of the order N c while the splittings are of the order 1/N C . One can easily argue that there 
are 0(1) negative corrections to M c \ which are in missing the present approach. Indeed, 
in some cases the so called Casimir energy was calculated for the solitonic solutions and 
was proven to be negative [25] . 

Hamiltonian Q has to be supplemented by a constraint which says that the allowed 
SU(3) representations 1Z = (p, q) have to contain states with hypercharge Y = 1 and the 
isospin of Y = 1 states is interpreted as (minus) baryon spin. The constraint Y = 1, which 
follows from the Wess-Zumino term, selects the representations of triality zero [7j . [o'] . [T2~] : 

8, 10, TO, 27, 35, 35, 64, . . . (5) 

for N c = 3. The success of the model is the prediction that the lowest baryonic states 
belong to the octet and decuplet representations of SU(3). The splitting between exotic 
(10, 27, etc.) and nonexotic representations depends on the strange moment of inertia I 2 , 
whereas the splitting between "standard" representations depends only on Ii. 

The wave function of baryon B = Y,I, J3 belonging to representation TZ = (p, q) and 
of spin S is defined as |2I 



\R,B,S) = VRi-f'D^ yrjjM) ( 6 ) 

where Y' = N c /3 and S = S, S 3 denotes spin and Q' is charge of the state (Y'JJ S ). 

are the SU(3) Wigner matrices in representation 1Z j22j. The quantization results in the 

relation: J = S and J3 = —S3. R denotes representation and its size as well: 

R = dim{n) = -(p+l)(q+l)(p + q + 2). (7) 
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Note that A is an SU(3) rotation matrix and therefore baryonic wave functions "live" 
in the SUfl avor (3)®SU S p in (3) space, where the second group is constrained to the SU(2) 
subgroup corresponding to spin. 

3. Symmetry breaking and the mass formulae 

Hamiltonian (JHJ) and the wave functions Q are identical in both Chiral Quark Model 
and Skyrme Model (with an obvious difference as far as model formulae for M c \ and are 
concerned). In the leading order in the symmetry breaking parameter m s (strange quark 
mass) and in the leading order of the N c expansion, the symmetry breaking Hamiltonian 
looks also identically in both models: 

H' = m cl - aD® (A) = m s a - m s oD { S (A) . (8) 

Indices a, b = 8 mean that Y = 0, 1 = 0, I3 = 0. Here a is related to the pion-nucleon 
sigma term: 

a = \^ (9) 
3 m 

m being the average mass of the up and down quarks. The classical mass m c \ = m s a is a 
constant which will be in what follows ignored, or more precisely included in the soliton 
mass M c i in Eq.©: M d -> M d + m ch 

(8) 

The symmetry breaking Hamiltonian D$ 8 is not reach enough to reproduce the mass 
splittings in the octet [22] • It has become evident that in order to fit the mass splittings 
one has go beyond the leading order either in N c or in m s . In 1984 Guadagnini [Z| 
introduced an additional term proportional to hypercharge, (3Y, whose coefficient (3 is of 
the order of 0(N®) and linear in m s . In the xQM still another nonleading term is present 
[26J, so that the perturbation Hamiltonian takes the following form: 

3 

H' xQM = -atDjg + (3Y + j=J2 D&Si (10) 

* i=l 

where a' = a + a, with a ~ 0(N C ) and a, [3, 7 ~ 0(N®) (of course all the constants in 
(fTUJ) are linear in m s ). Since the matrix elements of H' x q M between octet and decuplet 
wave functions © depend only upon two linear combinations of the three constants a', 
(3 and 7, the total number of parameters is 4; the two others being M d — and jj-, 
or equivalently M 8 and M w - average octet and decuplet masses. The above discussion 
shows that there is no way to extract I2 and, say, 7 from the spectra of the "standard" 
baryons. The quality of the fit to the mass spectra is very good and there exists a new 
(as compared to Gell-Mann-Okubo mass relations) sum rule known as the Guadagnini 
relation jjj: 

8(N + H*) + 3S = 11A + 8S*. (11) 

Another point of view was advocated by Yabu and Ando |28j who in 1988 proposed to 
diagonalize exactly breaking Hamiltonian (JBJ). By analytical methods developed in Ref. j 
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they reduced the problem to the harmonic oscillator with the frequency proportional to 
m s . They obtained a very good fit to the data, although their approach was criticized 
because at the level 0(m 2 s , m? s , . . .) there are in principle other possible symmetry breaking 
contributions to H' which they neglected. However, the same objection can be raised 
against H' q U and N c power counting, with an explicit example discussed for instance 
in Ref. |29j. Guided by the Yabu-Ando method we have observed in ^2] that already 
the second order perturbation theory in m s gives a good approximation to the fully 
diagonalized Hamiltonian. Moreover, as we will show below, the second order perturbation 
theory gives us a possibility to constrain the parameter I 2 responsible for the center of 
mass of the exotic antidecuplet. 

Let us observe that the symmetry breaking Hamiltonian 

H' SM = -aDlS(A) (12) 

does not change the quantum numbers of the state on which it acts; it can only change 
the representation: 

8«{ R) = (R',B,S\D®\R,B,S). (13) 
Therefore the first order correction to the baryon mass reads 

M bIr) = -«*£(*), (14) 
whereas the second order correction takes the following form 



R'^R 

where 



^R'R 



A R , R = C 2 (R')-C 2 (R). (16) 
Let us give for completeness the values A's: 

^To8 = 3, A 2 78 = 5, A2710 = 2, A 3510 = 6, A 8 jQ = —3, A 27 jq = 2, A 35 jq = 6, 

and the matrix elements 5 R , ^ : 



(17) 



5% = 


3 

10 ' 


= 


1 

20' 




i 2 = 


6 
100 


*S = - 


1 

10' 




1 

20' 




i 2 = 


4 

100 


°A — 


1 

10' 


K) 2 = 


0. 




i 2 = 


9 
100 


*i = - 


2 

10' 


(5l°) 2 = 


0. 


(si 7 : 


i 2 = 


6 
100 



for the octet, and for the decuplet: 



- h (^a) - m> (^a) - m> 

= 0; (^S*) = 128! (^| 5 *) = 896' 

Sg = -i = ± (*i) 2 = H, 

Sh° = 'I {^f = 0, = §:. 
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Finally for the antidecuplet we have: 



v e+; 


= o, 


(€v) 2 = 


72 

896 ' 


0$) = 


3 

640' 




81 

896' 


(^) 2 = 


8 

640' 


(^) 2 = 


72 
896' 


(€ /2 ) 2 = 


15 

640' 


m,f = 


45 
896 



JJ10 _ 2 /r8 \ 2 _ n 

°e+ — 8' \°e+) ~ u ' 

X10 _ 1 (XS \2 J_ 

U N_* 8' V U N*/ 20' \"JN~; 640' \"JN;/ 896' /iQ"i 

54 = 0, ^r = j- , ft**-* 
*%» = -i (^3 /2 ) 2 = ' 

Equations 0, (fT5j) and ()15jl form a complete mass formula for baryons in all allowed 
representations which depends upon 4 free parameters M ch Ii, J 2 and a, or equivalently 
on 

e = I 2 a 2 (20) 

a, Mg and Mio- By fixing these parameters from the nonexotic baryonic spectra we can 
predict masses in the exotic antidecuplet. For example the mass of B + reads: 

M e+ = M a+ 3 £- 2 - a -JL, (21) 

Similarly to the previous case ([TT]) the new sum rules can be derived by examining the 
null vectors of the m s correction matrix. For example the following mass relation 

132 AT 11„ 33^ 77 A 

N H H = — £ H A 22 

65 5 26 26 v ; 

is satisfied with a few promile accuracy. Similar sum rule holds for the decuplet: 

2627 A 213 , 1136 ,71 

A + S* = £* + —Q. 23 

594 22 99 27 v ; 

Somewhat awkward coefficients in Eqs. (|22|23|) are due to the SU(3) Clebsch-Gordan 
coefficients [27j entering formula ()15j) . 



4. Numerical estimates 

Question arises how to use the mass formula given as a sum of H of Eq.© M^ R) and 

MgL given by Eqs. (|14|15jl . In our early work from 1987 [12] we have estimated the value 
of parameter a and then fitted the remaining 3 parameters Ms,io and e by minimizing 
the functional 

F[M 8 M w a,e] = V (M$ [M 8 M 10 a, e] - Mf p ) 2 (24) 
18 B 

with a fixed. In 1987 the accepted value of the pion-nucleon sigma term was at the level of 
60 MeV [30J in contrast to the present value of 45 MeV [3T] used in Ref . jTSJ • The average 
light quark mass was taken as m = 5.5 MeV which gave a = 3.64. Taking m s ~ 200 MeV 
we got a ~ 720 MeV. By minimizing ([2*4^1 with respect to the remaining parameters we 



6 



Table 1: Parameters of the model and corresponding moments of inertia in MeV 



M 8 Mio a e l/I x l// 2 
fit I 1213 1503 720 1442 193 360 
fit 11 1230 1535 616 1824 205 284 



got the B + mass equal to 1534 MeV. This is fit no. I. In fit no. II we have allowed all 
parameters to be free with the result Mq+ = 1339 MeV. The parameters for the fits are 
given in Tabled 

Unfortunately the mass of + is quite sensitive to the choice of parameters. This 
is illustrated in Fig^a where we plot the results of the constrained fits obtained by 
minimizing ()24|) for fixed a. We see a rather steep rise of M e + with a. Two dashed 
vertical lines correspond to a of fits I and II and thin horizontal lines to the experimental 
masses. 




Figure 1: (a) Spectrum of the nonexotic baryons as a function of parameter a, together 
with a prediction for the mass of G + . Thin horizontal lines correspond to the experimental 
data, dashed vertical lines correspond to fit II (a = 615) and fit I (a = 720). (b) Spectrum 
of the exotic antidecuplet as a function of parameter a. 

Although there is rather large sensitivity of Mq+ to the choice of a, the conclusion 
one can draw from FigOa is that G + is relatively light. The full spectrum of the exotic 
antidecuplet as a function of parameter a is shown in Fig^b We see that below a = 600 
MeV the second order correction reshuffles the order of the antidecuplet states which 
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means that the perturbation theory becomes unreliable. As compared to the first order 
perturbation theory, the states are no longer equally spaced, with 9 + being significantly 
lighter than the other members of 10. It is interesting to note that if we fixed the value of 
I2 (or equivalently a) by the requirement that Mn* = 1710 MeV, we would get a = 738 
MeV corresponding to + mass of 1580 MeV. 

Finally let us present another estimate of I2 based on the model formulae for the 
parameters of Hamiltonian H in the first order of the perturbation theory. These pa- 
rameters depend on the soliton size, r , which enters the convenient variational Ansatz 
[32J of the soliton profile function P(r) (J2J): 



P(r) = 2 arctan 



rv 2 



(25) 



In practice one uses dimensionless variable Xq = eF^r®. Here e is a free parameter corre- 
sponding to the Skyrme term in the mesonic Lagrangian (JI} and F n = 186 MeV (for the 
details see e.g. Ref.[33]). In Ref.[33] we have calculated M cl , and Ji )2 in terms of x for 
the arctan Ansatz (|25|1: 

M cl = ^vr 2 ^(4* + -), (26) 
e 16 x 

I 1 = J F -n^(6xl + 25x ), (27) 

h = ^ F y^r(44 + 9x ). (28) 

By minimizing M c \ with respect to the soliton size we get we get Xq = a/15/4. It is easy 
to convince oneself that by fitting the octet-decuplet mass difference: 

M w - M 8 = = 231 MeV (29) 

we get e = 4.64. With this value of e we have 1/J 2 ~ 400 MeV and A% = 1762 MeV. 
Assuming that splittings in 10 are identical with the splittings in the ordinary decuplet 
(which is true in the first order of the perturbation theory for the breaking Hamiltonian 
flE} ) we get M©+ ~ 1460 MeV. 

5. Comparison with the Chiral Quark Model 

The authors of Ref. used the mass formula in the first order in m s with, however, 
nonleading terms in N c (see Eq. fjlOj) ). The reason was, as in our case, that the mass 
formula in the leading order in N c and in m s (JSJ) was unable to reproduce the nonexotic 
mass spectra with good enough accuracy. Systematic expansion of the rotational effective 
Hamiltonian in inverse powers of N c [2E] introduces 3 unknown free constants a', (3 and 
7 to H' of Eq.(|l()jl. However, only two linear combinations enter the nonexotic mass 
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splittings. Therefore the ordinary baryons allow to fix only 4 parameters M 8 , M 10 and 
the above mentioned two linear combinations of a', (3 and 7. In other words, in contrast 
to our approach, there is no way to fix I 2 without some further assumptions. Also the 
third linear combination of breaking parameters cannot be fixed. This is an important 
ingredient, because the exotic spectra are sensitive to other linear combinations of a', (3 
and 7 [13]. 

The main uncertainty due to I2 was removed in Ref. by the assumption that 
the "nucleonic" member of the antidecuplet was identified with the nucleon resonance 
N*(1710). The remaining freedom, illustrated in Fig|2la was removed by fixing m s = 150 
MeV and T* n N = 45 MeV. How big is the residual uncertainty due to the choice of m s and 
Etttv? We illustrate this in FigElb where we plot the spectrum of the antidecuplet states, 
for fixed N*(1710), as a function of a dimensionless parameter 

£71- AT 



s 



m. 



(30) 



150 50 

The choice of Ref. jT^j corresponds to S = 0.9. However, other choices of S are phe- 
nomenologically not excluded. For example in Refs.jHl] we have used parameterization 
with 7 = 0, as the explicit model calculations suggest which corresponds to S — 1.16. 
With this choice of S we get M®+ = 1560 MeV. This is very close to 1580 MeV obtained 
within the Skyrme model second order perturbation theory by fixing N* at 1710 MeV. 
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Figure 2: (a) Parameters a', (3 and 7 as functions of S constrained to fit the nonexotic 
baryon spectra, (b) Spectrum of antidecuplet in the Chiral Quark Model as a function of 
S. S = 0.9 corresponds to the choice of Ref. [13], whereas S = 1.16 corresponds to 7 = 0. 

Our original choice of m s and corresponds to S = 1.6, i.e. as seen from 

FigEla, to j3 = 0. Although in our approach 7 = as well, however, the second order 
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perturbation theory generates effects similar to 7 7^ 0. Figure Ela indicates that the good 
fit to the nonexotic spectra with S = 1.6 requires rather large a ~ 700 MeV in agreement 
with our fit I. 

6. Discussion and conclusions 

The above discussion shows that chiral models predict the existence of the low lying 
exotic baryons. The precise value of the lightest member of the exotic antidecuplet de- 
pends crucially on the estimate of the strange moment of inertia Ii- We have presented 
two kinds of estimates of 1 2 - The one based on the second order mass formula for baryons, 
allowed us to fix J 2 from the spectrum of the nonexotic states . We used the second order 
perturbation theory in order to fit the nonexotic spectra with high accuracy. With our 
original choice of parameters for the pion-nucleon sigma term E^jv — 60 MeV and the 
strange quark mass m s = 200 MeV we have predicted M@+ = 1534 MeV. This number is 
in striking agreement with the latest experimental evidence, however, one should not for- 
get that in fact the prediction of M e + depends rather strongly on the model parameters. 
If we used, as in Ref.jTH], the mass of N*(1710) to anchor the center of 10, we would get 
M®+ = 1560. A similar analysis has been recently done in Ref.[35 . 

The second estimate of I2 was based on the first order mass formula and on model 
expression for I 2 (|28p. With this method the mass of + comes out slightly lower; M@+ ~ 
1460. 

The existence of the light, exotic, positive parity baryons belonging to antidecuplet 
is a natural feature of the chiral models jS] , fOj [H] , [12] • From the beginning it was clear 
that the minimal quark content of such a multiplet must be q 4 q [O]. Quark models, how- 
ever, have difficulties in accommodating B + of positive parity (yet requiring experimental 
confirmation) unless strong interquark correlations are introduced jHHl > |HI1 • Perhaps the 
most striking difference between the quark and the soliton pictures consists in the fact 
that in the quark models 10 is inevitably accompanied by an exotic octet which does not 
appear in the soliton approach. With the discovery of + low energy QCD regained at- 
tention and the constituent quark vs. soliton interpretations of the light baryons compete 
in explaining the data. 

We would like to thank B. Jaffe, M. Polyakov, V. Petrov and K. Goeke for conversa- 
tions. We benefited from the e-mail exchanges with K. Hicks, D.I. Diakonov, M. Karliner 
and H.J. Lipkin. This work was partly supported by the Polish State Committee for 
Scientific Research under grant 2 P03B 043 24. 
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